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ABSTRACT 

Observations of the BO. 2 V star r Scorpii have revealed unusual stellar wind characteristics: 
red-shifted absorption in the far- ultraviolet OVI resonance doublet up to ~ +250 km s" 1 , and 
extremely hard X-ray emission implying gas at temperatures in excess of 10 7 K. We describe a 
phcnomenological model to explain these properties. We assume the wind of r Sco consists of 
two components: ambient gas in which denser clumps are embedded. The clumps are optically 
thick in the UV resonance lines primarily responsible for accelerating the ambient wind. The 
reduced acceleration causes the clumps to slow and even infall, all the while being confined 
by the ram pressure of the outflowing ambient wind. We calculate detailed trajectories of 
the clumps in the ambient stellar wind, accounting for a line radiation driving force and the 
momentum deposited by the ambient wind in the form of drag. We show these clumps will fall 
back towards the star with velocities of several hundred km s _1 for a broad range of initial 
conditions. The velocities of the clumps relative to the ambient stellar wind can approach 2000 
km s _1 , producing X-ray emitting plasmas with temperatures in excess of (1 — 6) x 10 7 K in bow 
shocks at their leading edge. The infalling material explains the peculiar red-shifted absorption 
wings seen in the OVI doublet. Of order 10 3 clumps with individual masses m c ~ 10 19 — 10 20 
g are needed to explain the observed X-ray luminosity and also to explain the strength of the 
OVI absorption lines. These values correspond to a mass loss rate in clumps of M c ~ 10~ 9 to 
10~ 8 M Q yr _1 , a small fraction of the total mass loss rate (M ~3x 10~ 8 M Q yr" 1 ). We discuss 
the position of r Sco in the HR diagram, concluding that r Sco is in a crucial position on the 
main sequence. Hotter stars near the spectral type of r Sco have too powerful winds for clumps 
to fall back to the stars, and cooler stars have too low mass loss rates to produce observable 
effects. The model developed here can be generally applied to line-driven outflows with clumps 
or density irregularities. 



Subject headings: stars: early-type — stars: individual (r Scorpii) — stars: winds — ultraviolet: 
stars — X-rays: stars 
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1. INTRODUCTION 

The star r Scorpii (HD 149438) is the MK standard for the BO V spectral type, although Walborn 
(1995, private communication) has classified it as BO. 2 V on the basis of its stellar wind lines. The star has 
served as a benchmark in studies of radiatively driven winds of hot stars. This is because it lies in a pivotal 
region of the HR diagram where stellar winds are making a transition from the massive and fast winds of 
the O stars to the winds of B stars that can just marginally be driven by line forces. 

In the Lamers & Rogerson (1978) study of Copernicus satellite spectra of early-type stars, r Sco was 
the latest spectral type to show the superionization stage of O VI. Stars with spectral types close to that 
of t Sco are also the latest to show the X-ray to bolometric luminosity relation L x /Lf, i = 10~ 7 that 
holds throughout the main sequence O spectral region. For stars of spectral type B1.5 V and later, this 
ratio of luminosities decreases by more than an order of magnitude, (Cassinelli et al. 1994). It is generally 
thought that the X-rays from hot stars arise from shocks embedded in the stellar winds (Owocki, Castor, & 
Rybicki 1988, hereafter OCR; Cooper 1994), so the decrease in the X-ray luminosity may be an additional 
consequence of the slower and less massive winds of near main sequence B stars. Massa, Savage, & Cassinelli 
(1984) found from a study of the UV lines that the wind properties of stars of spectral types BO V to B2 
V are particularly sensitive to the stellar composition, with winds ranging from non-detectable to quite 
strong for stars with increasing metal abundances. The stars near this spectral range are also commonly 
seen as emission line Be stars, for which the stellar rotation may produce wind compressed equatorial disks 
as suggested by Bjorkman & Cassinelli (1993) and Owocki, Cranmer, & Blondin (1994) or rotationally 
induced bistability as suggested by Lamers & Pauldrach (1991) and Lamers et al. (1999). 

While t Sco is a standard star, its stellar wind is anything but typical. The two most unusual 
properties of its wind, and those that we will focus on, are the peculiar red-shifted absorption seen in wind 
lines of high ionization species such as O VI, which are not present in the photosphere of this star, and its 
excessively-hard X-ray spectrum. 

The first of the peculiar wind properties of r Sco was discovered by Lamers & Rogerson (1978) in 
their analysis of Copernicus observations. They noted that the absorption produced by OVI and NV 
extended redward of line center by as much as +250 km s _1 . For stars with massive winds, the strong 
UV resonance lines have P Cygni-type profiles, which are noted for the absorption shortward of line 
center owing to scattering in the expanding wind. The absorption longward of the line center in r Sco is 
unusual. Absorption redward of line center could occur if the observed resonance line is also present in the 
photosphere of the star. However, the superionization stages of OVI, N V, and P V should not be present 
in the photosphere of r Sco, which has an effective temperature of <~ 31, 000 K. Two ways to explain the 
redshifted component are mentioned by Lamers and Rogerson (1978): turbulence in the flow near the 
base of the wind, and infall of highly-ionized material. Both phenomena could produce a positive velocity 
component in these ionic lines. 

The second peculiarity of r Sco is its very hard X-ray spectrum and large X-ray luminosity. In their 
R OS A T survey of near-main sequence B stars, Cassinelli et al. (1994) found that r Sco has a high X-ray to 
bolometric luminosity ratio: L x /L\y \ s=a 10~ 6 ; this is an order of magnitude larger than most O and early 
B stars. Furthermore, they found that the X-ray spectrum has an anomalously luminous hard component 
above 1 keV. The hard X-rays from this wind are unusual in that they require plasma temperatures 
in excess of 10 7 K. Cohen, Cassinelli, & Waldron (1997; hereafter CCW) have more recently presented 
observations of r Sco obtained with the ASCA satellite. These observations confirm the earlier ROSAT 
results requiring a very high temperature component to explain the observed X-ray spectrum. The ASCA 
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observations were fit with a multi-temperature equilibrium plasma model characterized by temperatures 
of T = (7, 12, and 27) x 10 6 K with emission measures of EM X = (3.5, 0.8, and > 0.3) x 10 54 cm" 3 , 
respectively. 

The very high temperature gas inferred from the X-ray spectrum of r Sco is not explained in the 
standard picture of X-ray production in early type stars (see discussion in CCW) . Numerical simulations 
of radiative shocks formed in the unstable stellar winds of OB stars suggest that these shocks are unable 
to produce temperatures much greater than a few times 10 6 K (OCR; Cooper 1994). The shock-jump 
velocities in these simulations are limited to roughly half the local stellar wind velocity. The presence of 
very hot gas in the wind of r Sco suggests that a different mechanism is responsible for at least some of 
the X-ray emission from this star. Though several groups have attempted to model the radiatively-driven 
outflow from r Sco, they were unable to reproduce the observed hardness of the X-ray spectrum (e.g., 
Cooper 1994; MacFarlane & Cassinelli 1989; Lucy & White 1980). 

In this paper we explore a phenomenological model to explain both the hard X-ray emission and the 
redshifted absorption of high ionization lines. Our model supposes density enhancements, or clumps, form 
in the smooth stellar wind and subsequently fall back towards the stellar surface, creating bow shocks at 
their leading edges. We investigate the trajectories of such clumps in a stellar wind and find infall to be 
achievable assuming the optical depths of these dense parcels of material are large enough to choke off the 
radiative line driving force. Using our calculated trajectories, we estimate the temperature of the shocked 
gas swept up by the clumps and comment on the ability of this model to provide the high-temperature 
gas implied by the observed X-ray emission (CCW; Cassinelli et al. 1994). While the numerical details 
presented here apply to the specific case of r Sco, the model developed in this paper is applicable to most 
clumped, radiatively-driven outflows. 

The details and assumptions of our model are discussed in §^|. In §|3| we derive the equation of motion 
for solitary clumps created within the wind and present the results of numerical integrations of clump 
trajectories. We describe the temperature spectrum of shocked plasmas in our model and discuss the 
observed X-ray emission in the context of this model in §|| In § 5] we briefly discuss why r Sco may be 
unusual among early- type stars, and we summarize our work in §5|. 

Table [l] lists the stellar parameters adopted throughout this work for r Sco. Empirical determinations 
of the mass loss rate of r Sco from UV lines gives discrepant results of 7 x 10~ 9 (Lamers & Rogerson 
1978), 7 x 10~ 8 (Gathier et al. 1981) and 1 x 10~ 9 M Q yr" 1 (Hamann 1981). These differences are due to 
differences in the adopted ionization fractions of the observed UV resonance lines of C III, C IV, N III, N V, 
O VI, Si IV, and P V. The terminal velocity of the wind of r Sco cannot be derived from the observations, 
because the UV resonance lines are unsaturated. Lamers and Rogerson (1978) derived a lower limit of 2000 
km s _1 . Because of these uncertainties we adopt the values predicted by the modified line driven wind 
theory (Kudritzki et al. 1989) with the force multiplier parameters taken from Abbott (1982): fc=0.156, 
a=0.609 and 5=0. 12. This gives a mass loss rate of 3.1 x 10~ 8 M Q yr _1 and Voo — 2400 km s _1 , in 
agreement with the range of empirically derived parameters. 



2. THE MODEL 

The winds of early type stars are driven by resonance line scattering of radiation in the stellar envelope. 
This radiative driving mechanism is inherently unstable (Lucy and Solomon 1970; Lucy 1982; OCR). Small 
perturbations on an otherwise monotonic wind velocity structure tend to grow and steepen into shocks. 
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These instabilities have been well studied using one-dimensional hydrodynamical models (see, e.g., Cooper 
1994, Feldmeier 1995), which show that the resulting wind is a mixture of rarified hot material and colder, 
dense regions. These studies also suggest the shock velocity jumps have maximal amplitudes of order one 
half the local wind velocity. The relationship between the velocity amplitude of a shock, Av, and the 
temperature behind the shock, T s , is given by 



In the case of most main sequence O stars and OB supergiants, the shock jump velocities are indeed strong 
enough to produce gas at the temperatures suggested by the observed X-ray emission (T 3 x 10 6 K; 
Chlebowski, Harnden, & Sciortino 1989; Cassinelli et al. 1994). The X-ray spectrum of r Sco, however, 
implies gas with temperatures in excess of 10 7 K (CCW). Numerical models of B star winds cannot produce 
shocks with such high temperatures. Velocity jumps of ~ 300 km s" 1 are typical in numerical simulations 
(Cooper 1994), while velocity jumps of more than 700 to 800 km s _1 are required to reproduce the hard 
X-ray emission from t Sco. 

While the one-dimensional shock production problem has been well studied, little attention has yet 
been given to the production of shocks and the structure of the wind in two or three dimensions. We assume 
that shocks and density enhancements in the stellar wind are relatively localized phenomena whose sizes 
are dictated by the local sound speed and instabilities that may tend to break up large sheets, such as the 
Rayleigh- Taylor instability. An individual shock is assumed to cover a small solid angle when viewed from 
the star. The one-dimensional simulations assume spherical symmetry; our assumption is that the spherical 
shocks in these simulations would fragment relatively easily. Cassinelli and Swank (1983) have suggested 
such a scenario to explain the lack of variability observed in their Einstein Solid State Spectrometer 
observations of OB supergiants. A stellar wind containing an ensemble of N shocks that continuously 
form and dissipate can produce a non-detectable variability of order \/N /N. The time variability of single 
spherical shocks in the wind of r Sco is discussed by MacFarlane & Cassinelli (1989). 

A schematic of our model is shown in Figure [l]. The basic picture of the stellar wind that we adopt 
is one in which discrete density enhancements, or clumps, are embedded in an otherwise smooth stellar 
wind. These dense clumps of material have shock interfaces with the faster moving ambient wind. The 
most extreme of these clumps are dense enough to be optically thick to radiation in the strongest resonance 
lines, though still optically thin to UV continuum radiation. In addition to their treatment of stellar 
winds with organized velocity gradients, the theory of radiative acceleration on material with no internal 
velocity gradient was developed by Castor, Abbott, & Klein (1975, hereafter CAK), and we use their 
parameterization of the line-driving force. The equations governing this acceleration are given in f^. 
Assuming a CAK parameterization of the line driving force, the net radiative acceleration on optically thick 
parcels of material with no internal velocity gradient can be relatively small (see §[!]) . Gravity dominates the 
motion of these optically thick clumps, which fall inward, creating bow-shocks at their leading edges as they 
plow supersonically through the outflowing wind. The velocity structure of the ambient wind is assumed to 
be relatively unaffected by such clumps. They do not shadow a significant amount of the overlying wind, 
since they are optically thin to continuum radiation and the line opacity is significantly redshifted relative 
to the opacity driving the unperturbed wind. 

As the clumps fall inward, they shock the out-flowing gas of the ambient (inter-clump) wind. The 
temperature of gas the immediate post-shock region is proportional to the square of the shock jump velocity 
as in Equation ([!]). Since the relevant velocity for the calculation of the shock temperatures in these 
bow shocks is the relative velocity between the clumps and the ambient wind, the temperatures needed 
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to produce the hard X-ray component observed by ROSAT and ASCA are readily obtained. The X-ray 
producing shocks adjacent to the leading edge of the infalling clumps produce very highly-ionized material, 
thus accounting for the red-shifted high ions seen in absorption (O VI, NV, etc.). 

A related model was considered by Lucy & White (1980, hereafter LW), in the first paper attempting 
to explain the source of hot star X-rays as shock-heated gas. They suggested that clumps are radiatively 
driven outward while the rest of the wind is pushed along by the drag of these objects. In our model the 
clumps move outward (or inward) with a velocity less than that of the ambient wind, so the shocks occur 
on the inner face of the clump. This placement of the shocks is more consistent with the results of OCR, 
who find that shocks tend to be inward facing with low-density, high-speed wind material being driven into 
slower and denser shocked regions. However, we retain the clump aspect of the LW model. 

Let us assume the clumps are self-contained parcels of gas confined by the ram pressure of the wind 
through which they fall. Following LW, we assume a balance of dynamical pressures and write the average 
density of a clump, (p c ), falling through the wind as 

(p c )=p w (r) {ai+ f pUj2) . (2) 

a c 

In this expression p w is the wind density, a function of the distance r from the star, a w and a c are the 
thermal velocities in the wind and clump, respectively, and to is the relative velocity between the wind and 
the clump. We have ignored the effects of turbulence within the clumps but have included the constant 
C p = 0.29 as a scale factor to account for the range of angles at which the wind hits the bow-shock (see 
LW). We can see from this that once a clump begins to decelerate, thereby increasing oj, the oncoming 
wind confines it more effectively making it denser as it falls towards the star. This further decreases the 
ability of the radiation field to drive the clump outwards. If the mass of material accumulated into a clump 
is specified, then the volume can be readily obtained from the density derived using Equation (0). 

This model has the desirable attribute of producing both the red-shifted high-ion absorption and the 
anomolously hard X-rays. Several independent lines of evidence suggest that clumping can be important 
in the winds of Wolf-Rayet (e.g., Lepine, Moffat, & Henriksen 1996) and OB supergiant winds (e.g., Sako 
et al. 1999; Eversberg, Lepine, & Moffat 1998). It is within reason that near main sequence early type 
stars may also contain localized regions of significantly higher than average densities. The exact mechanism 
for the formation of these dense clumps is not rigorously approached in this work, but it is worthwhile to 
explore the consequences of such a model before a more thorough hydrodynamical study is undertaken. 



3. INFALLING MATERIAL 

The main driving force for winds from massive stars is the radiative acceleration associated with 
atomic absorption line opacity. The most important of the driving lines are the ultraviolet resonance lines 
of abundant ionic species. Thus, the driving force is not only dependent upon the elemental abundances in 
a wind, but also upon the wind's ionization balance (MacFarlane, Cohen, & Wang 1994; Abbott 1982). The 
equation of motion describing the trajectory of a parcel of gas with respect to the star is 



dv c _ GM. 
~~dt ~ 



, (l-T e -T L )+f d , (3) 



where v c is the velocity of the clump of material, M* the mass of the star, r the distance of a clump from its 
center, and fd is the force per unit mass the ambient wind exerts upon the clump in the form of a "drag." 
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The terms T e and T l are the ratios of the radiative to gravitational acceleration for electron scattering and 
resonance line scattering, respectively. The drag force, fd, is zero if there is no relative motion between a 
clump and the ambient wind. In this treatment we have not allowed for a thermal pressure gradient in the 
wind, which is negligible in the supersonic portion of such an outflow. 

For the ambient stellar wind, fd = 0. The requirement for the material to be accelerated away from the 
star is dv/dt > 0, implying T e + Tl > 1, i.e., that the radiative acceleration is greater than the gravitational 
acceleration. A dense clump of material formed in the wind that has T e < 1 may decelerate and 

possibly fall back to the surface of the star, plowing its way through the out-flowing wind as it falls. 
However, even for clumps with T e + Tl < 1, the drag force, fd, may be great enough to push the clump 
outwards in the stellar wind. To determine if infall is physically reasonable, we shall develop more fully 
each of the terms in Equation (0) . 



3.1. Radiative Acceleration 

The radiative driving force on a parcel of wind material is considered to have two components: radiative 
driving via electron scattering and resonance line opacities. The radiative acceleration due to electron 
scattering of the stellar continuum is written 

9e - ^; (4) 
c 

here a e is the electron scattering opacity, which is a e = axh/ l^eiTiH ~ 0.34 cm 2 g" 1 for r Sco, F is the 
frequency- integrated stellar flux at the location of the material, and c is the speed of light. 

A proper treatment of the radiative driving of stellar winds owing to absorption line opacity has been 
given by CAK (see also Abbott 1982). Their theory paramaterizes the line driving force on a wind having 
an organized velocity gradient dv/dr. In CAK theory, the radiative acceleration due to line scattering, g^, 
is written 

gL = ^l M (t)=gfM(t). (5) 
c 

where g r e ei is for a reference electron scattering opacity of 0.325 cm 2 g _1 . So g^ is simply the reference 
electron scattering acceleration modified by a multiplicative factor. This factor, -M(i), the so-called "force 
multiplier," is a function of t, an optical depth parameter that is independent of line strength. The force 
multiplier is defined as 

M(t) = kt~ a (jf) , (6) 

where N e is the electron density in units of 10 11 cm~ 3 , W is the dilution factor and k, a and S are the 
CAK parameters. These parameters have been determined from studies employing extensive atomic line 
lists and depend upon the temperature of the absorbing material (see, e.g., Abbott 1982). However, the 
driving force for the density enhancements we will consider here is slightly different from the typical CAK 
driving force for the ambient wind. 

In the CAK treatment of expanding atmospheres or winds, t is the electron scattering optical depth over 
the distance in which the wind velocity increases by the thermal velocity. Hence t is inversely proportional 
to the local velocity gradient for outflows with monotonically increasing velocity laws (CAK). In the case of 
parcels of gas with no organized velocity gradient, such as those being considered here, the optical depth 
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parameter t is the integral of the electron scattering opacity over the path length / through the region or 
clump [Equation (5) of CAK]: 

t = a e I pdr. (7) 



Approximating the integral as (p c }Vc^ 3 , where V c = m c /p c is the volume of the clump, Equation (Q) may 
be rewritten, using Equation as 

,2\ 2/3 



t » ffem yV(r) a/3 fi + ^- 



(8) 



assuming a c ~ o ffi . This latter assumption has the desirable consequence that the density of a clump will 
approach that of the ambient wind as the relative velocity, uj, approaches zero. 

The ratio of the radiative to gravitational accelerations in CAK formalism is written 

given L, the luminosity of the star. We will assume the values of the CAK parameters tabulated by Abbott 
(1982) for stellar winds at 30,000 K (k — 0.156, a — 0.609 and 5 = 0.12) are appropriate for use in the 
clumps considered here. These are the same as those used for estimating the mass loss rate. Given the 
high densities predicted for the clumps and the presence of a hard X-ray source at their leading edges, the 
ionization balance in the clumps could be different than that of the ambient wind. Thus while the CAK 
parameters appropriate for the clump may be different than those used here, we adopt these parameters in 
order to proceed. In the end, the choice of CAK parameters will likely not matter too much so long as the 
optical depth through a clump is very high. Equation (^) gives an upper limit to the radiative acceleration 
because it assumes that the flux from the star reaching a parcel of gas is unattenuated by the layers closer 
to the star. For a clump moving outwards at lower velocities than the ambient wind, the photons producing 
the line radiation pressure on the dense clump may already have been absorbed or scattered by gas lower 
in the wind with the same outward velocity as the clump. We will ignore this effect here, which is expected 
to be small. 



3.2. Drag force 

The standard drag force per unit mass on a clump of area A c moving through an ambient medium 
with a relative velocity lu is given by 

f d = lc DPw (r)u 2 ( — ). (10) 
2 \m c ) 

In this treatment m c is the mass of the clump, and p w {r) is the density of the ambient medium (in this 
case the ambient stellar wind). The drag coefficient, Cd, dictates the efficiency of momentum transfer to 
the clump, and the factor of x j i is conventionally included. This expression describes the energy density 
impinging on the clump integrated over the area of the clump. For most of this work we will assume 
Cd = 2, i.e., perfect transfer of momentum from the ambient wind to the clump. This is the most stringent 
assumption for determining if infall is viable. 

2/3 

The cross section of the clump may be approximated A c « V c , giving 

\C D p w uj 2 



(m c (p c )*) 1/3 ' 
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Using Equation (g), and assuming a c rs a w and ui 3> a c , we rewrite this as 

''-M^r) ■ (12) 

The density of the wind at a distance r from the surface of the star is calculated assuming mass conservation: 

M 

Pw[r) = A^My (13) 

where M is the mass loss rate (see Table |l|) . The velocity of the wind v w (r) is calculated using a (3 velocity 
law with (3 W = 0.8 (Groenewegen & Lamers 1989): 

v w (r)=v 00 (l-^j . (14) 

with the values of Voo and i?* given in Table @. 

The magnitude of the drag force can be relatively small. We estimate its magnitude using Equations 
(H) and |FJ) with the data given m Table |: 

A -66 oar- ( ] L)" /3 ( T ^-)" , \. W ./ 3 , („, 

where we have assumed to — v w , which in this equation is in units km s _1 . If we define A = fd/g, where g 
is the force of gravity per unit mass, we can rewrite Equation (^|) as dv c /dt = —g(l — T e — Tl — A). We find 

A , 6xl0 -3(l_) 4/ *(_-_)- I/a „„ (r) V», m 

where v w is again in km s _1 . We can see that very near the star, the drag force is relatively small compared 
with the force of gravity because the velocity between the clump and the ambient wind is smaller than at 
larger radii. This prediction is based on the assumption lo m v Wl but for clumps producing temperatures in 
excess of 10 7 K, our numerical calculations discussed in § |j.3| suggest the drag force is much more significant 
(relative to gravity) at larger radii. For clumps with masses logm c ~ 20.0 (in grams) the outward drag 
force is comparable to the force of gravity at large distances (r > 10i?*) from the stellar surface. For lower 
mass clumps, however, the drag force can play an important role in the dynamics at intermediate distances 
from the star (r ~ 5R*). 



3.3. Numerical Results 

We have derived the individual components of Equation (||), and can hence obtain trajectories in 
position-velocity space for clumps in the wind of r Sco. Equation (|^) shows the dynamics of such clumps 
depend on the relative velocity, to, between the clump and the ambient wind: lo = v w (r) — v c , where v w (r) 
and v c are the wind and clump velocities, respectively. We have numerically integrated the equation of 
motion for parcels of gas with various initial conditions since it does not have a analytically tractable 
solution. The initial conditions specified in our approach are the initial distance r Q from the center of the 
star, the initial relative velocity u> , and the mass m c of the clump (we assume mass conservation). By 
specifiying lu and r Q , we are also then specifying the initial velocity of the clump, given Equation (|l4|). We 
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assume the value of lo is characteristic of the jump velocities typically seen in numerical simulations, i.e., 
Av <J v w /2 (Cooper 1994). The dependence on m c comes in through Equations (||) and (|l^). the former 
affecting A4(t) in the radiative acceleration. 

Our numerical integrations of the equation of motion show two types of clump trajectories: outflow 
and infall trajectories. We term those trajectories in which the clumps are driven from the star, though 
they may initially decelerate, outflow trajectories. Though these clumps are driven outward, they always 
have velocities less than the terminal velocity of the wind. Those that decelerate and eventually descend 
towards the star follow infall trajectories. The latter type is the most interesting, given the redshifted high 
ion absorption line observations. It is not clear, however, that outflow and infall trajectories are mutually 
exclusive. There may exist a mix of trajectories based upon a distribution of initial conditions for the 
clumps. Figure displays the velocity of representative outflow and infall trajectories as a function of 
distance from the star. Also shown is the adopted (3 velocity law of the wind. 

To demonstrate the sensitivity of these trajectories to initial conditions, we have compiled a series 
of trajectories that are different in only one of the three possible input conditions. Figure |^ shows the 
trajectories of clumps with r = 1.20 i?» and tu = v w /2 — 331 km s" 1 ; the masses of the clumps in this 
figure are logm c = 18.0, 19.0, 19.5, and 20.0 (m c in grams). The low-mass clumps are driven outward 
through the wind while the high-mass objects relatively quickly stop and start to descend towards the 
star. The low-mass clumps continue outwards because of their lower column density (and correspondingly 
low tl), allowing for more effective radiation driving. Lower- mass clumps are also more effectively pushed 
outwards by the drag force provided by the on-rushing stellar wind, though this is usually minor compared 
with the force provided by resonance line scattering of radiation. Figure ^ shows trajectories for clumps 
having the same mass and initial velocity relative to the wind [(logm c , oj ) = (19.0, 250 km s -1 )] but with 
starting positions r a = 1.10, 1.15, 1.20, and 1.25 and in Figure ||| we show trajectories for clumps having 
a mass logm c = 19.0 starting at r a = 1.20 i?» with initial velocities relative to the wind of uj = 100, 200, 
300 and 400 km s^ 1 . Clumps with larger ui tend to fall in more readily than those with low tu , as parcels 
of gas with smaller initial distances from the stellar surface, r Q , are also more likely to fall inwards. We 
note that the factor C p — 0.29 included in Equation (Q) is not well known. If this constant were larger, the 
clumps in our models would be more effectively confined and have higher densities. This would make it 
easier for a clump of given initial conditions to fall inwards. 

In Table ^| we give several properties for trajectories with a range of initial conditions. Included in the 
table are the initial conditions w G , r D , and logm c . Also given is the lifetime T of the clump, in hours, the 
final velocity of the clump in our calculations, the maximum relative velocity between the clump and the 
stellar wind, uj m ax, and the corresponding maximum shock temperature. For the outflow trajectories, the 
final velocity is the velocity of the clump at a distance of 10i?» from the surface of the star, and the lifetime 
is the time it takes to reach this point. 

It is clear from Figures that conditions for infall can be met for a broad range of initial conditions. 
Our hypothesis that the infall of material may be responsible for the observed red-shifted component of 
the high-ionization lines is thus reasonable. Table || shows that many combinations of initial conditions 
can yield infall trajectories. Infall velocities of ~ 300 km s _1 (as observed in the red wings of the O VI 
lines) are found for clumps with masses in excess of ~ 10 19 grams. Clumps originating at larger distances 
from the star reach higher infall velocities. Clumps following infall trajectories spend more time near their 
"turning point," i.e., where v c ~ 0, than close to the photosphere where their infall velocity reaches its 
maximum. The observed O VI profiles require this, showing absorption consistent with a Gaussian optical 
depth distribution centered at v c ~ km s _1 relative to the photosphere. 
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Typical lifetimes for clumps following infall trajectories are ~ 3 — 10 hours, while outflow trajectories 
take 15 — 30 hours to reach r ~ 10/2* . In general clumps with initial distances r Q ^ 1.2/?* with 
uj = 0.5v w (r o ) require masses of order logm c > 19.1 to achieve infall. The highest-temperature gas is 
produced for clumps that reach the largest distances from the stellar surface. Thus lower mass clumps 
starting at relatively large distances from the star produce the hardest X-rays. One can see from Table g 
that the trajectories followed by high-density clumps with logm c = 19.5 tend to produce slightly harder 
X-rays than those with logm c = 18.0 or 20.0. We believe the most likely clump masses for r Sco are in the 
range logm c = 19.0 to 19.5 (in grams) given their ability to infall and produce very hot gas. 

For infall trajectories u max occurs near the turn-around point in a cloud's trajectory. Specifically, the 
largest velocity jumps are usually found immediately after a clump of material has begun to fall inwards. 
At this point the parcel of gas is at the maximum height above the star, and the ambient wind velocity 
is larger than at any other point in the clump's trajectory. Though clumps may accelerate as they near 
the star, the decrease in the wind velocity with distance from the stellar surface offsets this increase in the 
infall velocity of the material. 



4. X-RAY EMISSION 

4.1. Temperature Distribution 

Having shown that the interpretation of the red-shifted high-ion absorption as infall is a physically 
reasonable one, we now approach the anomolous X-ray spectrum for r Sco. We have speculated that 
the observed X-rays are a by-product of the interaction of infalling clumps with the outgoing wind. The 
temperature T s for gas in a section of a bow shock can be calculated by substituting v s cos 6 for Av 
in Equation ([l]), where v s is the shock-jump velocity (here the relative velocity between the wind and 
a clump, w) and 9 is the angle at which the material enters the shock, relative to a ray normal to the 
shock's surface. The temperature of gas is a function of position along the shock's leading edge, via the 
9 dependence. Calculating the true run of temperatures requires knowledge of the exact shape of the 
bow-shock. This is an inherently difficult problem and is beyond the scope of this work. When discussing 
the temperatures produced by our infalling parcels of gas, we will typically give results that assume 9 = 0°, 
and, following LW, values corresponding to the shock at 9 — 30°. Truth may be somewhere between these 
values for the majority of the shocked gas. 



The trajectories calculated in § |3.3| permit us to estimate the bow-shock temperature at each point of 
a clump's history given its velocity relative to the ambient wind. It is clear that the infall velocities of the 
trajectories displayed in Figures |^-^ can exceed ~ 300 km s _1 , sometimes approaching ~ 500 km s _1 near 
the base of the wind. As mentioned in § |3.3| , the relative velocity between the wind and the infalling clumps 
can be much higher than these values. The maximum relative velocities in the trajectories for clouds having 
logm c = 18.0, 18.5, 19.0, and 19.5 in Figure | are u; max « 1460, 1660, 1295, and 1202 km s" 1 . These 
velocities imply maximum shock temperatures of order T s sa (31, 40, 24, and 21) x 10 6 K, respectively, 
well in excess of the temperatures required to match the CCW and Cassinelli et al. (1994) fits. Even the 
outflow trajectories provide very large velocity jumps, and hence large shock-jump temperatures. If wc 
assume Av = u> cos 30° we find T s « (5 — 10) x 10 6 K. This is still quite close to the temperatures required 
to explain the ASCA observations of CCW. 

Figures ^| and ^ show two trajectories as well as the model temperature distribution of the 
shock-heated gas. The initial conditions for the trajectory shown in Figure || are (logm c , r a , to ) = 
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(19.0, 1.3 i?*, 371 km s" 1 ), while those for Figure | are (logm C! r a , uj ) = (19.0, 1.5 R*, 498 km s" 1 ); 
in each case lo = v w /2. Changing the initial height above the stellar surface between these two figures 
causes the clump to flow outwards in the wind. The maximum shock temperatures expected from these 
two trajectories are near 5 x 10 6 K. The insets show the expected relative temperature distribution over 
the lifetime of the clump. These distributions have been weighted by p w (r) 2 for each point in the clumps' 
trajectories to approximately account for the expected differences in the X-ray emission measure at different 
points in the wind. We show both the distribution assuming 8 — 0° (solid histogram) and 9 — 30° (dotted 
histogram) in these figures. 

The temperature distributions shown in Figures |5| and |TJ show that the trajectories of dense parcels 
of gas derived in §|| provide large enough shock-jump velocities to produce gas in excess of 10 7 K. This 
treatment is very approximate. For instance we have neglected the cooling of the gas in the bow shocks 
around the clumps. The gas that is heated to T ~ 10 7 will eventually cool through the lower X-ray emitting 
temperatures, providing for enhanced emission in the soft X-ray bands over that shown in Figures ^ and 
||. Though the temperature distributions peak at the same temperatures, the temperature distribution for 
the infall trajectory (Figure ||) is narrower than that of the outflow trajectory (Figure [3]), which shows a 
strong tail towards higher temperatures. The ratio of the emission measures of X-ray emitting gas observed 
by CCW at (1.2 and 2.4) x 10 7 K is <J 3.5. Both the outflow and infall trajectories have similar displayed 
in these figures produce X-ray emitting material at (1.2 and 2.4) x 10 7 K in a ratio of 1.9:1. This value 
is consistent with the observational constraints, though the effects of radiative and adiabatic cooling will 
affect the final distribution of X-ray emission. 

In all likelihood a population of clumps in the wind of r Sco will have a range of initial conditions, 
providing a range of X-ray producing temperatures and X-ray luminosities per clump. The distribution 
of intial conditions, the true shape of the bow shock at the leading edges of the infalling clumps, and the 
cooling of the post-shock gas will all determine the true temperature and luminosity distribution of X-ray 
emitting gas. The derivation of these factors is beyond the scope of this paper. Our goal is to show that 
the production of very hot gas, in excess of 10 7 K, is possible in the wind of r Sco through the infall of 
dense material. Our dynamical modelling of the dense clumps in §|^ has shown that infall can occur. The 
discussion presented in this section shows that the infall can create the requisite conditions for producing 
this very hot gas 

4.2. X-ray Luminosity and the Number of Clumps 

The X-ray luminosity predicted in our model is also reasonable. The mechanical luminosity of the wind 
L w ~ VzMw^ is a factor of ~ 780 greater than Lx for this star. Even if we assume the appropriate velocity 
for use is not the full terminal velocity, but rather ~ 1500 — 2000 km s _1 , similar to the observed wind 
velocities at distances typical of our model clumps, we still have L w ~ 300 — 600Lx- Thus the ensemble 
of clumps need only convert < 1% of the wind mechanical luminosity to X-ray luminosity to match the 
observed Lx- 

We estimate the available luminosity, L Cl for X-ray production from a single clump as L c ~ 1 /2/? tu w 3 j4 c . 
The trajectory in Table || with initial conditions (logm c , r a , oj ) — (19.0, 1.20 i?*, 286 km s _1 ) reaches 
a maximum height of ~ 1.5 i?» with a relative velocity u> ~ 1000 km s _1 at this height. This clump will 
have an area A c ~ 3 x 10 20 cm 2 at its maximum height, and we find L c ~ 5 x 10 29 ergs s _1 . Thus of 
order ~ 150 individual clumps, assuming perfect efficiency, are required to produce the observed X-ray 
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luminosity of the star, Lx ~ 7.3 x 10 31 ergs s^ 1 (Cohen, Cassinelli, & MacFarlane 1997; hereafter CCM). If 
we assume an efficiency of ~ 10 — 20% for conversion of kinetic energy to X-ray luminosity in the ROSAT 
bandpass (see, e.g., Wilson & Raymond 1999), then the total number of clumps required by the X-ray 
observations is of order 10 3 . The observed luminosity of the hard component is significantly lower than the 
total {Lx ~ 1.8 x 10 31 ergs s _1 ; CCW), in qualitative agreement with the temperature distributions shown 
in the insets of Figures || and ^|. 

We can estimate a lower limit to the number of clumps from the lack of (short-term) variability in the 
observed X-rays (CCM). The population of clumps that exist in the wind at any one time (N) should be 
large enough to produce X-ray variations less than yf~N /N <; 0.04, or N > 600 clumps. Assuming 600 clumps 
of mass logm c ~ 19.0 with lifetimes similar to those in Table ^, T ~ 5 — 11 hours, this yields a creation rate 
of mass in clumps of M c ~ (2 — 5) x 10 -9 M Q yr _1 . Assuming half of the clumps actually escape in the 
wind (follow outflow trajectories), the mass loss rate in clumps is M c ~ (1 — 3) x 10~ 9 M© yr _1 , or ~ 10% 
of M. 

Another constraint on the number of clumps residing in the stellar wind of r Sco comes from the O VI 
absorption line profile. The A1038 A profile shows ~ 10% absorption near v ~ km s^ 1 . Assuming the 
O VI optical depth for an individual clump is much larger than unity and, as above, that the clumps are 
at a distance r ~ 1.5 i?* with mass log?7j c f=a 19.0, approximately N c ~ 1900 clumps are required to cover 
10% of the total solid angle of the star. This corresponds to a mass loss rate of M = (0.8 — 2.0) x 10~ 8 M© 
yr _1 for assumed lifetimes T ~ 5 — 11 hours. If the efficiency of mechanical to X-ray luminosity from the 
previous paragraph is ~ 10%, these two estimates of the number of clumps are in agreement. It is worth 
noting that the infall trajectories are naturally weighted to provide more absorption near v c ~ km s _1 , 
i.e., near the upper-most point in the trajectory, than at larger infall velocities. This can be seen in Figure 
U where we have plotted a dot every 1.25 hours along the clump's trajectory. The high density of dots near 
the turning points implies that a random sample of clumps drawn from from this trajectory will be heavily 
weighted towards low velocities, in agreement with the Gaussian shape of the O VI absorption. Detailed 
line-profile calculations will be the subject of a future work. 

Our calculations suggest that the required number of clumps in the wind of r Sco is of order 10 3 , with 
M c approximately a few x 10~ 9 M Q yr _1 , consistent with the adopted mass loss rate of this star. These 
rough numbers can vary significantly depending on the assumptions made, however. At the present time 
there is no clear evidence for an inconsistency with the observational constraints. Further work will need to 
place more stringent constraints on the allowable parameters of this model. 

5. WHY r SCO? 

Though our simple model can explain several of the observational aspects of the stellar wind of r Sco, 
there are several questions and details that we have not approached. Principal among these is what makes 
t Sco so special? Why do most hot stars not show the observed peculiarities of this star? This is may 
be due to the special location of r Sco on the main sequence where the mass loss rate steeply drops with 
T e ff. Main sequence stars hotter than B0, i.e., the main sequence O stars, have mass loss rates considerably 
higher than r Sco. For instance, the mass loss rate of 09 V stars with T c g ~ 35, 000 K, L ~ 10 5 Lq, and 
M ~ 25M Q are of order 10 -7 M Q yr _1 (predicted by the CAK theory), compared with less than 10~ 9 or 
10~ 10 M Q yr _1 for B1.5 and later stars. This is observed in the sudden disappearance of the P Cygni 
profiles of the UV resonance lines in stars between 09.5V and B0.5V (see the atlas of UV line profiles of 
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Snow et al. 1994). The mass loss rates of B1.5 V stars are typically a factor of ~ 10 below r Sco (e.g., the 
stars in CCM), with B2 star mass loss rates often being a factor of 100 lower than that adopted for r Sco. 

We speculate that the prevailing conditions for stars near the position of r Sco in the HR diagram 
are just right not only for the formation of dense clumps of material within the wind, but also for their 
ability to infall. Stellar winds near the spectral type of r Sco are only marginally able to be driven by the 
radiation of the underlying stars. Clumps with no organized velocity gradient formed in such winds are less 
likely to be pushed outward than in the winds of earlier-type stars. 

To investigate the differences between r Sco and earlier-type O stars, we have calculated trajectories 
for high-density clumps in the wind of the 09 V star 10 Lacertae. We adopt stellar parameters given in 
Brandt et al. (1998) and Snow et al. (1994) in our calculations: = 10 5 L Q , M„ = 25 M , R* = 8.3 R Q , 
M = 1.6 x 10~ 7 M yr" 1 , T eff = 36,000 K, and Voo = 2.6v esc = 2520 km s" 1 . We note that like r Sco, 
10 Lac has a low projected rotation velocity (vsinz ss 40 km s^ 1 ). Wc find no infall trajectories for initial 
distances from the stellar surface r a ;> 1.03 i?*, initial relative velocities u> <; 0.5v w (r o ), and masses in 
the range logm c < 21.0. This is primarily a function of the increased luminosity to mass ratio of the star, 
though the extra drag force plays a limited role. This explains the lack of red-shifted absorption components 
of O VI in the spectrum of 10 Lac. 

Our calculations show that clumps in the wind of 10 Lac will have shock jump velocities in excess of 
1000 km s _1 , producing temperatures in excess of 10 7 K. This is very close to the value of 0.5v w found for 
shocks in "normal" stellar winds. The emission from (outflowing) clumps in the wind of 10 Lac will be 
indistinguishable from the normal X-ray flux due to shocks in unstable stellar winds of O-stars. 

The lack of redshifted absorption profiles in the UV lines of main sequence stars of spectral type B0. 5 
and later is due to drastic decrease in mass loss. A simple scaling of the optical depth of the O VI lines 
by a factor 0.1 compared to r Sco, due to a ten times lower mass loss rate, makes these lines undetectable 
in the Copernicus spectra, though few such stars were observed with Copernicus. The significantly lower 
mass loss rates of later-type stars may not provide the requisite conditions for formation of an ensemble 
of clumps like those thought to exist in the wind of r Sco. In general we would expect fewer, lower mass, 
clumps to form in the winds of later type stars. Such clumps would be less dense than their counterparts in 
the wind of r Sco, hence more likely to flow outward. In the winds of later-type B stars, with significantly 
lower densities and terminal velocities, the clumps would produce lower temperature and lower luminosity 
X-rays, perhaps undetectable or indistinguishable from the expected shock emission. 

6. SUMMARY 

We have outlined a model to explain the two peculiar observational properties of the stellar wind of 
r Sco: infall of high-density clumps surrounded by very highly-ionized gas that gives rise to the red-shifted 
absorption wings of the O VI lines, and strong shocks that produce the hard X-ray spectrum. 

Our model hypothesizes that density enhancements in the stellar wind lead to the formation of dense 
clumps of material, confined by ram pressure, that decouple from the dynamical flow of the ambient stellar 
wind. These clumps have very large optical depths in the UV resonance lines responsible for driving the 
stellar wind. The resulting radiative acceleration of the clumps is insufficient to overcome the gravity, so 
they slow down and fall back towards the star. As they do so, their interaction with the on-rushing stellar 
wind produces a bow shock, with shock temperatures reaching several times 10 7 K. 
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We have presented detailed dynamical models of density enhancements in the wind of r Sco, which 
show that infall can be achieved for a wide range of initial conditions. Using these trajectories we have 
shown that the shock-jump velocities can approach 2000 km s _1 . The anomolously hard X-rays observed 
from t Sco by Cassinelli et al. (1994) and CCW are produced in these bow shocks at temperatures of 
~ (1 — 5) x 10 7 K. The infalling high-ionization material observed in the OVI doublet by Lamers & 
Rogerson (1978) is material in and around the infalling clumps that is ionized by the hard X-rays produced 
in the bow shocks at the clumps' leading edges. 

Our calculations show that clumps with masses in the range of 10 19 to 10 20 g will fall back if they 
are formed at radii <; 1.3 i?* (see Table 2). Only the clumps formed above 1.1 i?* produce shocks with 
temperatures in excess of 10 7 K. Clumps formed this high have maximum infall velocities of order 500 
km s _1 . This is higher than the observed maximum redshift of 250 km s -1 , but the clumps spend very 
little time at such high velocities. Even clumps with a large final infall velocity will give the maximum 
contribution to the the UV line profile at small velocities, because they spend much more time near their 
turning point, where v ~ 0, than at high velocity just before they disappear into the photosphere. This 
explains why the OVI absorption reaches a maximum near km s _1 . 

We estimate that ;> 10 3 clumps of mass 10 19 — 10 20 g are needed to explain the observed 
OVI absorption lines and the hard X-ray flux. This corresponds to a clump mass loss rate of 
^clumps ~ l-5x 10~ 9 M yr _1 , or ~ 0.1 M win d- In making these estimates we have assumed several 
efficiency factors whose values are not well known at present and would require detailed modelling to be 
improved. For instance, (a) in comparing the predicted X-ray luminosities with observations we have 
assumed a perfect conversion of ram pressure heating to X-ray emission; (b) we have neglected the radiative 
cooling of shocks which may result in a reduction of the predicted X-ray flux from an individual clump; 
and (c) we have assumed that the clumps are optically thick in the O VI transitions, and that the intrinsic 
line widths of the clumps is sufficient to account for the breadth of the observed absorption. A reduction in 
these "efficiencies" may result in an increase in the required number of clumps. This will be the subject of 
a subsequent study. 

We explain the special situation of r Sco in terms of its location on the main sequence where the mass 
loss rate drops by about factor 100 from 09 V to Bl V. This can be tested and our model can be improved 
by high S/N observations of the far UV resonance lines of OVI, S VI, and P V of a number of 09.5 V to 
Bl V stars with the Far Ultraviolet Spectroscopic Explorer (FUSE). This satellite will cover the wavelength 
range 905 — 1197 A and will observe a large number of early- type stars. Our model would suggest that 
stars with spectral types similar to r Sco should also show infall in the O VI doublet at AA1032 and 1038 A, 
which will be well observed with FUSE. 

The dynamical model developed in this work can be generally applied to clumps in line-driven outflows. 
Though it seems that the properties of stars near r Sco on the HR diagram may be the best for producing 
infalling clumps with recognizable observational consequences, clumping may occur in many types of 
early- type stellar winds, and may also be important in outflows from extragalactic sources (e.g., AGNs and 
QSOs). The dynamical model described here can be applied to all of these cases (as discussed for the case 
of 10 Lac above). 

We thank D. Cohen and A. Cole for many conversations regarding this project. JCH recognizes support 
from a NASA Graduate Student Researcher Fellowship under grant number NGT-5-50121. JPC recognizes 
support from NASA under grant NAG5-2854. 
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Tablc 1. Stellar Parameters for r Scorpii 



Quantity 


Value 


Source 


Spectral Type a 


BO. 2 V 


Walborn(1994) 


Distance (pc) 


132 


Hipparcos (Perryman et al. 1998) 


V (mag) 


2.82 


Hipparcos (Perryman et al. 1998) 


Nhi.ism (cm" 2 ) 


2.9 x 10 20 


Diplas & Savage (1994) 


T eff (K) 


31,400 


Kilian (1992) 


R* (-R©) 


6.5 


Snow et al. (1994) 


M. (M Q ) 


15 


Snow et al. (1994) 


w esc (km s" 1 ) 


917 


Snow et al. (1994) 


L» (ergs s _1 ) 


1.4 x 10 38 


Snow et al. (1994) 


L x b (ergs s" 1 ) 


7.3 x 10 31 


Cohen et al. (1997) 


vs'mi (km s" 1 ) 


20 


Uesugi and Fukuda (1982) 


M c (Mq yr- 1 ) 


3.1 x 10~ 8 


This work 


Woo (km s _1 ) 


2400 


This work 



a The Bright Star Catalog (Hoffleit & Jaschek 1982) has the 
classification B0 V. 

b We have adjusted the value given in Cohen et al. (1997) to reflect 
the revised distance used here. 

c Theoretical calculations based on modified CAK theory, with 
values taken from Abbott (1982) and the using the fitting formula of 
Kudritzki et al. (1989) 
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Table 2. Properties of Clump Trajectories 



log m c 


To 


LOo 


v w (r ) 




Ub, final 


, i d 


rp e 

J. s max 


[g] 

LOJ 


[R*\ 

L ^ J 


[km s _1 l 

L J 


[km s _1 l 

L J 


[hrs.l 

L J 


[km s _1 l 


[km s _1 l 

L J 


[10 6 K] 

L J 


18.0 


1.10 


176 


352 


6.3 


-230 


730 


8 


18.0 


1.25 


331 


662 


22.4 


+790 


1460 


31 


19.0 


1.10 


176 


352 


2.8 


-270 


580 


5 


19.0 


1.15 


235 


470 


4.5 


-340 


820 


10 


19.0 


1.20 


286 


572 


6.8 


-400 


1055 


16 


19.0 


1.25 


331 


662 


11.2 


-450 


1295 


24 


19.0 


1.30 


371 


742 


25.9 


-490 


1590 


36 


19.0 


1.50 


498 


996 


25.1 


+540 


1670 


40 


19.5 


1.10 


176 


352 


2.5 


-280 


560 


5 


19.5 


1.15 


235 


470 


3.8 


-355 


785 


9 


19.5 


1.20 


286 


572 


5.5 


-415 


1000 


14 


19.5 


1.25 


331 


662 


7.9 


-470 


1202 


21 


19.5 


1.30 


371 


742 


11.8 


-515 


1405 


28 


19.5 


1.50 


498 


996 


33.7 


+385 


1820 


48 


20.0 


1.10 


176 


352 


2.3 


-285 


555 


4 


20.0 


1.20 


286 


572 


4.8 


-430 


970 


13 


20.0 


1.30 


371 


742 


9.3 


-530 


1340 


26 


20.0 


1.50 


498 


996 


58.5 


+300 


1990 


57 



a Velocity of the stellar wind at the initial radius, r Q , of the clump. 

b Lifetime of the clump in hours. For outflow trajectories, this is the time it takes 
the clump to reach r = 10i?». 

c Final velocity of the clump. Trajectories with Vb, final < are inflow trajectories, 
while final > are outflow trajectories. For outflow velocities Vb, final is the velocity 
at r = lOi?*. 

d Maximum relative velocity between the wind and the clump. 
°Maximum shock temperature obtainable given the value of uj ma x- 
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A MODEL FOR 
i SCORPII 



2. Bow skock forms 
as amkient wind 
flows past klok. i 



1. Blok forms from 
amkient wind. ^ 



-t;i 
c 




3. T < 1; Blok slows 
and stops. 








^.-ray emission 
(10 7 Kgas) 



Infalling high -ions 
(OV1, NV, PV) 




4. Blok falls 
kack onto 
tke star. 



Fig. 1. — Illustration of the production of X-rays and infalling matter in the form of clumps in the stellar 
wind of r Sco. The clumps are formed from the ambient stellar wind. Because the radiation force is lower 
than that of gravity (r e + Tl < 1), the clumps slow, stop, and then fall back onto the star. The hard 
X-rays are formed in a bow shock where the wind collides with an infalling clump. The X-rays produce 
the superionization seen in the UV spectra by way of the Auger ionization mechanism. The wind is also 
assumed to have the shocks typically asssociated with hot stars; these produce the soft component of the 
X-ray spectrum of r Sco. 
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Fig. 2. — Trajectories of clumps forming in the wind of r Sco having different initial conditions. The velocity 
scale is relative to the surface of the star, and the distance from the stellar surface is given in units of stellar 
radii. Also shown is the (3 velocity law of the ambient stellar wind (solid line). The starting distance from 
the star and velocity relative to the ambient wind for these clumps are (r , u> ) = (1.25 R*, 331 km s _1 ). 
The four different trajectories represent clumps formed with masses logm c = 18.0, 18.5, 19.0, and 19.5 (m c 
in grams). The masses are labelled for each trajectory. Note the low-mass clumps are pushed outwards 
through the wind, following outflow trajectories. The higher mass clumps are able to fall back onto the star 
and follow infall trajectories. 
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Fig. 3. — As Figure but for clumps of constant mass and initial velocity relative to the wind [(log m c , u ) = 
(19.0, 250 km s -1 )] with varying initial radii [r = 1.10, 1.20, 1.25, and 1.30 J?*]. The initial radii are 
labelled for each trajectory. 
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Fig. 4. — As Figure [| but for clumps of constant mass and initial radii [(logm c , r a ) = (19.0, 1.25 i?*)] 
with varying initial relative velocities [ui — 150, 250, 300, and 400 km s -1 )]. At this position in the wind, 
1 /2V W — 331 km s _1 . The initial relative velocities are labelled for each trajectory. 



- 23 - 




Fig. 5. — The main plot shows a trajectory for a clump with the initial conditions (logm c , r a , lu ) = 
(19.0, 1.30 371 km s _1 ). This clump follows an infall trajectory. The trajectory is marked with a dot 
every 1.25 hours. Also shown is the wind velocity calculated with a (3 velocity law (dot-dashed line). The 
inset shows the relative distribution of temperatures in the shock-heated gas formed as the clump follows this 
trajectory. This distribution is weighted by p w (r) 2 for each point along the trajectory to roughly account 
for the dependence of the X-ray emission on density. The shock-jump velocities, Aw, are calculated from 
the relative velocity between wind and the clump, u), such that Av — lucosQ. The solid distribution of 
temperatures assumes 9 — 0° , while the dotted distribution uses 6 = 30°. 
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Fig. 6. — As for Figure ||, but for the initial conditions (logm c , r , u ) = (19.0, 1.50 i?*, 498 km s 1 ). 
This clump follows an outflow trajectory. 



